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Long range electrostatic, induction and dispersion coefficients including terms of order R~ 8 have 
been calculated by the sum over states method using time dependent density functional theory. 
We also computed electrostatic moments and static polarizabilities of the individual diatoms up 
to the octopole order using coupled cluster and density functional theory. The laboratory-frame 
transformed electrostatic moments and van der Waals coefficients corresponding to the alignment 
of the diatomic molecules were found. We use this transformation to obtain the coupling induced 
by an external DC electric field, and present values for all XY combinations of like polar alkali 
diatomic molecules with atoms from Li to Cs. Analytic solutions to the dressed-state laboratory- 
frame electrostatic moments and long range intermolecular potentials are also given for the DC 
low-field limit. 



I. INTRODUCTION 

Advances in the formation of ultracold absolute ground 
state polar alkali diatoms [l], Q open up avenues into 
many branches of the physical sciences. For chemical 
physics, applications of polar diatoms range from preci- 
sion spectroscopy 0,13, to the study 0-Q and control Q 
of cold chemical reactions. Other areas of physics benefit 
from the use of polar molecules, such as condensed mat- 
ter physics [Tfj| , with the search for novel quantum gases 
fill and phases [l^ . Furthermore, dipolar gases have 
been the subject of much interest from the quantum in- 
formation community [l3l - [T5j . and ideas of atom optics 
(e.g. using evanecent wave mirrors [IH) have been gen- 
eralized to polar molecules [13, EH ■ The recent achieve- 
ments in molecular alignment and control flfjl |20| may 
also allow to take advantage of the unique properties and 
possible control provided by ultracold polar molecules. In 
addition, there is growing interest in reactions of alkali 
diatoms to form tetramer structures [2ll - |23j with reason- 
able dipole moments and rich molecular structures, which 
could offer good candidates for quantum computing with 
dipoles (24[. In each of these applications it is crucial 
to accurately describe the intcr-molccular interactions, 
themselves dominated by their long range behavior [25[ 
at the low temperatures found in these systems. Be- 
cause of the weakness of the long range intermolecular 
forces as compared to the chemical bond, and the range 
of nuclear coordinates and phase space involved, it is ad- 
vantageous to consider alternative methods of modeling 
the intermolecular potential other than ab initio quan- 
tum chemical calculations. 

A standard approach to describing the long range in- 
teraction potential between two molecules, in the limit 
that the wavefunction overlap between the molecules is 
negligible, is to expand the interaction energy into three 
distinct components, 

Eint = E c \ + Ei n d + -Edisp- (1) 

Here E e \, E^ and -Edisp are the permanent electrostatic, 
induction (permanent-induced electrostatic) and disper- 
sion energies. Each of these terms can be perturbatively 



expanded in an asymptotic van der Waals series, 

E hn = J2 C " R ~ n - ( 2 ) 

n 

The coefficients C„ are in general angular dependent, and 
can be computed in several ways. In this work we expand 
the intermolecular electronic interaction operator in a 
multipole expansion (2(| , and then use first- and second- 
order perturbation theory to calculate the van der Waals 
coefficients. Several papers have discussed the isotropic 
i?~ 6 interactions of homonuclear alkali diatoms using 
both the London approximation [27j and time depen- 
dent density functional theory (TD-DFT) dHH]. The 
isotropic and anisotropic contributions have been inves- 



FIG. 1. (a) Schematic representation of an aligned diatomic 
molecule. Classically, the molecule precesses on a cone of an- 
gle about the electric field F, with {cos 9) describing the av- 
erage orientation of the molecule: its dipole moment T> points 
towards F for (cosO) > 0, and in the opposite direction for 
(cos 8) > 0. The alignment, (cos 2 8), describes the tightness 
of the rotational cone, (b) Lab-fixed frame molecular interac- 
tion geometry in the presence of an external field, where Of is 
the angle between the field and the vector R joining the two 
molecules. 
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tigated using configuration interaction [3(| and TD-DFT 
[311 ] to compute van der Waals coefficients through R~ s . 
However, systematic research on the heteronuclear alkali 
diatoms is limited to the R~ 6 isotropic van der Waals 
coefficients for the LiX (X=Na,K,Rb,Cs) species Q. To 
date, the only heteronuclear anisotropic van der Waals 
coefficients available in the literature are for KRb and 
RbCs [32| and limited to R~~ 6 dispersion forces. In this 
paper we present a systematic study of the isotropic and 
anisotropic van der Waals interactions through order R~ s 
of the heteronuclear alkali rigid- rotor diatoms in their ab- 
solute ground state as a continuation of our work on the 
homonuclear species [3l|. Also included is the transfor- 
mation of the long range interaction potential from the 
molecule-fixed (MF) frame to the laboratory-fixed (LF) 
frame for use in molecular alignment computations. Af- 
ter a brief description of dressed state diatomic molecules 
in Sec. UH we review the sum over states method of cal- 
culating van der Waals coefficients in Sec. IIII1 In Sec. 
IIVI the transformation to, and matrix elements of, the 
lab-fixed frame van der Waals interaction potential are 
described. Analytic expressions of the low-field field cou- 
pled electrostatic moments and van der Waals coefficients 
are also provided. The ah initio methodology is outlined 
in Sec. [V] and we conclude in Sec. [VI] with a discussion 
of our numerical results. 



II. DRESSED STATE DIATOMIC MOLECULES 

The orientation and alignment ( (cos 9) and (cos 2 9) [37| 
respectively as illustrated in Fig. [TJa)) of polar molecules 
can be achieved through several mechanisms, the most 
direct of which is the coupling of rotational states by a 
polarizing external DC electric field, F. Increasing the 
strength, F, of the external electric field increases the 
number of rotational states coupled, thus tightening the 
orientation of the molecule in a cone of angle 9 about the 
orientation of the field. To account for this rotational 
coupling adiabatically, we expand the dressed state rota- 
tional wave function of the molecule as a superposition 
of field-free symmetric top states 

| JMQ) = all I J M ty , (3) 

J,M 

where the symmetric top states are given in term of 
Wigncr rotation matrices D J _ M _ n (a, (3, 7) [38| 

\JMSl) = (-l) M - n (^^) Dt M _ n (a,P, 7 ), (4) 

where (a, /?, 7) arc the Euler angles of the molecule, and J 
the total angular momentum quantum number with pro- 
jections M in the laboratory frame and ft onto the molec- 
ular axis. The expansion coefficients a J M dictate the levels 
of mixing between the different rotational states, and can 
be solved for by diagonalizing (J' M'fl'\H\JMn) . Here 



H is the symmetric top and dipolc-ficld Hamiltonian 

H = B(J 2 + Jl)-VFcos9, (5) 

where J is the angular momentum operator, J z is the 
angular momentum projection on the z axis, B is the 
molecular rotational constant, T> is the dipole moment 
of the molecule, and 9 is the angle between the exter- 
nal electric field of magnitude F and the molecular axis. 
The coefficients a J M (F) then depend on the strength F of 



TABLE I. Center of mass multipole electrostatic moments, 
{Qto), of all the ground state heteronuclear alkali diatoms 
through cesium evaluated at the equilibrium bond lengtlQ r e . 
The variable R q denotes the distance where the R~ 5 electro- 
static term overcomes the dipole-dipole R~ 3 contribution. All 
values are presented in atomic units. 



System Method 


r a 


(Q10) 


(Q20) 


<Q30) 


Rq 


LiNa 


CCSD(T) 
CCSDT. 


5.45 
5.43 
5.45 


0.20 
0.22 
0.21 


10.07 


-47.33 


95 


LiK 


CCSD(T) 
CCSDT. 


6.27 
6.21 
6.27 


1.39 
1.39 
1.38 


6.07 


-59.99 


15 


LiRb 


CCSD(T) 

VCI^ 

CCSDT. 


6.50 
6.48 

a. ten 


1.63 
1.63 


2.76 


-62.41 


13 


LiCs 


CCSD(T) 
CCSDT. 


6.93 
6.82 
6.93 


2.15 
2.17 
2.11 


-2.29 


-49.88 


10 


NaK 


CCSD(T) 

ccsd(t;£ 
vc£ 


6.61 

6.592 

6.49 


1.12 

1.156 

1.09 


10.56 
10.60 


-26.54 


19 


NaRb 


CCSD(T) 
VCI^ 


6.88 
6.84 


1.35 
1.30 


6.94 


-56.00 


16 


NaCs 


CCSD(T) 

vc£ 


7.27 
7.20 


1.85 
1.83 


2.49 


-60.45 


12 


KRb 


CCSD(T) 

vc£ 

reP~ 


7.69 
7.64 
7.7 


0.25 
0.23 
0.30 


15.14 


-69.09 


109 


KCs 


CCSD(T) 
VCI^ 


8.10 
8.02 


0.75 
0.76 


13.00 


-105.70 


38 


RbCs 


CCSD(T) 
VCI b 


8.37 
8.30 


0.49 
0.40 


15.88 


-50.28 


60 



a r e values are taken from experimental results where available, 

see Deiglmayr et al. [33j and references therein. 
b Ref. [3. 
c Ref. |T 
d Ref. 231 . 

e Ref. 35| performed a four component Dirac-Fock valence bond 
calculation in calculating the dipole moment. 
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TABLE II. Multipole static polarizabilities, a u > m , and isotropic van der Waals dispersion coefficients, W^qqq , up to order 
n — 8 of all the ground state alkali diatoms through cesium evaluated at the equilibrium bond lengths r e listed in Table |T] All 
values are presented in atomic units. 



System 


Method 


a 

Qll(l_ 






^220 


Q221 


Q222 


"6000 


"sooo 


LiNa 


PBEO 

VCP 

CCSDT^ 


300.0 
347.6 


185.5 
181.8 


223.7 
237.0 
237.8 


9418.9 


7035.6 


3356.2 


3.279[3] 
3.673[3£ 


4.982[5] 


LiK 


PBEO 

vcp_ 

CCSDT^ 


455.1 
489.7 


261.8 
236.2 


326.3 
320.7 
324.9 


24164.4 


15899.8 


5939.6 


5.982[3] 
6.269 [3£ 


1.378[6] 


LiRb 


PBEO 
VClf, 
CCSDT^ 


445.5 
524.3 


256.1 
246.5 


319.2 
339.1 
346.2 


27815.3 


18110.7 


6359.2 


6.193[3] 
6. 323 [3^ 


1.583 [6] 


LiCs 


PBEO 
CCSDT^ 


525.2 
597.0 


289.1 
262.5 


367.8 
374.0 
386.7 


38723.9 


24996.3 


7935.8 


7.700[3] 
7.712[3]1 


2.297[6] 


NaK 


PBEO 
VCP_ 

CCSD(T^ 


472.7 
529.2 


280.6 
262.3 


344.6 
351.3 
363.8 


16572.0 


13035.0 


6739.5 


6.818[3] 
6.493[3£ 


1.268 [6] 


NaRb 


PBEO 
VGlf, 


504.6 
572.0 


285.3 
280.3 


358.4 
377.5 


25217.0 


17771.7 


7547.5 


7.688[3] 


1.790 [6] 


NaCs 


PBEO 
vrT d 


587.3 


323.2 


411.2 


37633.3 


25245.7 


9444.6 


9.453[3] 


2.641 [6] 


KRb 


PBEO 
VCI?. 


729.6 
748.7 


420.9 
382.9 


523.8 
504.8 


36974.1 


27588.8 


13100.9 


1.349[4] 


3.385[6] 


KCs 


PBEO 
VCI?. 


836.7 
822.3 


468.6 
425.62 


591.3 
571.1 


56372.9 


38791.7 


16262.9 


1.657[4] 


5. 038 [6] 


RbCs 


PBEO 
VCI d 


901.0 
904.0 


502.0 
492.3 


635.0 
602.8 


48325.3 


36401.8 


18619.8 


1.884[4] 


5. 188 [6] 



a Note that the parallel and perpendicular static dipole polarizabilities, ay and ctj 
b a = h(ctu + 2ckx) is the average static dipole polarizability. 



correspond to il'm = 110 and 111 respectively. 



Ref. 
d Ref. 
e Ref. 



Ref. 



s Ref. 



evaluated using CCSD and the Tang-Slater-Kirkwood formula [36 
2711 evaluated using the London formula. 



the field. While theoretically simple, this process can be- 
come experimentally challenging for molecules with small 
dipole moments or rotational constants due to the large 
external fields required for strong alignment. 



state is of interest. As such we present our alignment in 
terms of an applied external static field and, where prac- 
tical, the number of strongly coupled rotational states. 



An alternative to simply increasing the static field 
magnitude is to add a separate polarizing laser field [39[ 
that directly couples the rotational states of the molecule. 
However, to achieve both alignment and orientation con- 
trol, time-dependent nonadiabatic effects are introduced 
into the dressed state wavefunction |2(} ■ For the purposes 
of this work the investigation and inclusion of these nona- 
diabatic effects arc unimportant as only the final dressed 



III. ANISOTROPIC LONG RANGE 
INTERACTIONS 

Given a linear molecule in the Born-Oppenheimer ap- 
proximation (no nuclear motion), at any given configura- 
tion the orientation of each molecule can be described by 
the vector f, = fa) with the relative position between 
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TABLE III. LiX (X=Na,K,Rb,Cs) calculated CCSD(T) elec- 
trostatic and TD-DFT dispersion+induction van der Waals 



( 1 2\ 

coefficients, W^£ L M , for unique combinations of L1L2M. 
All values are presented in atomic units and calculated at the 
equilibrium bond length r e listed in Table [I] and [n] denotes 
xHT. 



n L\ L 2 M 


LiNa 


LiK 


LiRb 


LiCs 


Electrostatic: 


W w 

yv nL 1 L 2 M 








3110 


-7.076[0] 


-3.799[0] 


-5.170[0] 


-9.094[0] 


3111 


0.038[0] 


1.900[0] 


2.585[0] 


4.547[0] 


4210 


-5.893[0] 


-2.547[1] 


-1.431 [1] 


1.324[1] 


4211 


1.965 [0] 


8.489[0] 


4.771 [0] 


-4.412[0] 


5220 


6.086 [2] 


2.276 [2] 


5.284[1] 


2.568[1] 


5221 


-1.353[2] 


-5.058[1] 


-1.174[1] 


-5.707[0] 


5222 


8.453[0] 


3.161[0] 


0.734[0] 


0.357[0] 


5310 


3.674[1] 


3.306[2] 


2.887[2] 


3.942[2] 


5311 


-9.186[0] 


-8.266[1] 


-7.218[1] 


-9.855[1] 


Dispersion+induction: W 


(2) 

iL\L 2 M 






6000 


3.289[3] 


7.243[3] 


7.254[3] 


1.062 [4] 


6200 


4.036[2] 


0.4 f » |ZJ 




7QQM 1 


6220 


1.768 [2] 


1.180 [3] 


1.094 [3] 


2.567[3] 


6221 


-3.929[1] 


-2.621[2] 


-2.430[2] 


-5.705[2] 


6222 


4.911[0] 


3.277[1] 


3.038[1] 


7.131[1] 


7100 


1.075[3] 


-4.898[2] 


-5.906 [3] 


-3.740[4] 


7210 


1.460[2] 


-7.641 [3] 


-8.641 [3] 


-3.110[4] 


7211 


-2.433[1] 


1.273[3] 


1.440 [3] 


5.184[3] 


7300 


7.411[2] 


4.567[3] 


1.057 [3] 


-9.020[3] 


7320 


4.259[2] 


-7.013[2] 


-2.682[3] 


-1.687 [4] 


7321 


-7.099[1] 




4.4 / U[ZJ 


Z.o±Z[oJ 


7322 


5.070[0] 


-8.349[0] 


-3.193[1] 


-2.008[2] 


8000 


5.586[5] 


1.539 [6] 


1.715[6] 


2. 722 [6] 


8200 


3.552[5] 


1.270 [6] 


1.534 [6] 


2. 793 [6] 


8220 


9.460[4] 


5.702[5] 


6.695[5] 


1.636[6] 


8221 


-1.406[4] 


-8.078[4] 


-9.379[4] 


-2.272[5] 


8222 


1.234[3] 


4.607[3] 


4.642[3] 


9.978[3] 


8400 


3. 896 [4] 


4.574[4] 


6. 320 [4] 


1.371 [5] 


8420 


2. 706 [4] 


7.284[4] 


8. 377 [4] 


2.485[5] 


8421 


-3.666[3] 


-9.822[3] 


-1.124[4] 


-3. 320 [4] 


8422 


2.001 [2] 


5.287[2] 


5. 965 [2] 


1.740 [3] 



the molecular center of mass denned as R = (R,6,(p). 
Here 9i is the projection angle of fj on R, <f>i is the 
projection angle of fj on the x axis and (R,9,(j)) are 
the spherical vector components of R. Due to the ro- 
tational invariance of the interaction energy between two 
molecules, it can be separated into a series of radial and 
angular basis functions 

£int(ri,f 2 ,R) = E LlLaL (R)A LlLaL (fi,f2,R)- (6) 

Li ,L 2 ,L 

Here El 1 l 2 l(R) are purely radial functions for a rigid- 
rotor and A[ yl L 2 [ y (fi 1 r2 1 R) is an angular basis which, 
when R is oriented along the z axis, can be expressed as 



A Ll L 2 L{fl,f 2 ,R) 



min(L-L,L 2 ) 

e vl i 1 l 2 l 

M=0 
M /„„„ n \ r>M 1 



x P£(cos0i)P£(cos0 2 )cos[M(0i - 2 )], (7) 



where 



n M 

r IL 1 L 2 L 



(-1) M (2 - 6m,o)(L 1 M;L 2 - M\L0) 



(Li — M)!(L 2 — M)\ 
(Li + M)!(L 2 + M)\ 



1/2 



(8) 



{L\M\L2 — M\L0) is a Clebsch- Gordon coefficient, and 
Pjf (cos 0) is an associated Lcgcndrc polynomial. The 
radial functions El 1 l 2 l(R) can be evaluated using first- 
and second-order perturbation theory by expanding in 
terms of the electronic multipole operators Qem = 
Y^i z i r iCe m {fi), where the sum is over all charges, Z{ is 
the charge at each z'th center, r\ is the distance from 
each i'th charge to the center of mass and Ci m (fi) is a 
Racah spherical harmonic [38j . 

Following the standard approach (2(| |U |4l|, the 
first- and second-order interaction energy for two linear 
molecules can be expressed as 



Eint(R, 9l,02 



E 

n,Li,L 2 ,M 
jM / „ a \ oM 



1?" 



where 



x (cos 9i)P% (cos 62) cos[Aty], (9) 



i) 



nh\L 2 M 



(-1)^(2-<W )- 



(L 1 +L 2 y. 



{L x + M)\{L 2 + M)\ 
x (Oi|Q Ll o|Oi)(0 2 |Q L20 |0 2 ) (10) 



is the first-order electrostatic contribution, where |0j) is 
the electronic ground state of molecule i. In Eq.©, 



W. 



(2) 



- w {2 ^ IS) 



' nL 1 L 2 M ~ " ' nh x L 2 M yy nL 1 L 2 M 

contains the second-order contributions from dispersion, 



m/ (2,dis) / m _ Me'^e' 

W nL 1 L 2 M\ K ) - ^LiL 2 M ^i+^i 



I2A 



E 



',+i 2 +i' 2 +2,n 



rpOiki rp0 2 k 2 



-, (12) 



and induction, 

m/ (2.ind) , m _ fia'^t^* 



\+e 2 +e'+2,n 



t 1 ,i 1 
I2A 



k 2 ^0 



rp0 2 k 2 



£ fc 2 - £o 2 



+ (1^2) . (13) 
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The scalar coupling coefficient C 



LxL 2 M 



2 is given WCl as 



ClIm 2 = (-l)^((2Li + 1)!(2L 2 + l)!) 1 ^ 
"(24 + 24 + l)!(2ii + 2^ + 1)!] 1/2 ^ M 



(24)!(2f 1 )!(24)!(24; 



x (£ 1 +£ 2 0;£' 1 +£' 2 0\L0) { £ 2 



i\ Li 
£> 2 L 2 

£i+£ 2 £[+£ 2 L 



(14) 



the symbol between curly brackets being a Wigner 9-j 
symbol [IH, and T^/ L is the coupled monomer multi- 
pole transition moment defined as 

Tf$ Ll = Y^( i\Qeirn\ki){k i \Q ei - m \Oi)(l i m-,£' i - m|40), 

m 

(15) 

where the indices fci go over ground and excited states 
of the i'th molecules electronic wavefunction \ki) [4^|. It 
is convenient, when discussing molecular properties, to 
work with the uncoupled dynamic multipolc polarizabil- 
ity: 



a«'m(w) = ^ 
fe#o 



(ffc ~ £o)(0i\Qe m \ki)(h\Qt'- m \0i) 
(efc - f-o) 2 - 



(16) 



The zero frequency limit of Eq. (1161) represents the static 
multipole polarizability. 



IV. DRESSED-STATE VAN DER WAALS 
INTERACTION 

A. General Expressions 

To consider the interactions between rigid-rotor linear 
molecules dressed by an external electric field it is neces- 
sary to first transform the van der Waals interaction en- 
ergy from the molecule-fixed frame (MF) to the lab-fixed 
frame (LF). The lab-fixed frame van der Waals interac- 
tion can be generally expressed by referring to Eq.© 
and removing the constraint on Eq.© which specified 
that R is aligned to the z-axis. The angular basis can 
then generally be expressed [4l[ as 



E 

mij ,m L2 ,m L 



f L t L 2 L 
\m Ll m L . 2 m L 



x Y LimLi (h)Y L2mL2 (f 2 )Y LmL (R), (17) 



where Yf mf (r) is a spherical harmonic and (:::) is a 
Wigner 3-j symbol [38| ■ Because of the change in angu- 

lar basis, it is necessary to rccouplc the radial W^£ L M 
functions. This can be done readily by integrating Eq.(j9|) 
over the angular phase space: 



1 



/> 7T />7T 

' d6t / d6»2sin(6»i)sin( 
o Jo 



where A 1 , 



Sir Jo 

x A' LiL2L {e 1 e 2( f>)E^l{R,e 1 ,e2,cp), (is) 

h,02,4 > ) is Eq.lfTTf projected onto the 



molecule-fixed frame and is given by [42 



A' 



min(Li,X/2) 



/ 2L+1 
V 2tt 



1/2 



x J2 (-l) m (2-5 m , ) 



rn=0 



L\ L 2 L 

m —m 



x e Ll m(0i)e Lam (02)cos[m0]. (19) 



where 0z m (£?) are normalized associated Legendre poly- 
nomials. The resulting integrand has the solution 



M 



I Li L 2 L 

M -M 



(Li + M)\(L 2 + M)l\ 1/2 
(Li - M)\{L 2 - M)\ 



2L + 1 



(2Li + l)(2L 2 + l) 



1/2 



E 



w 



(2) 



m) (20) 



R r ' 



The adiabatic dressed state basis for two molecules at 
large separation is given in terms of the product of each 
molecule dressed rotational wave functions 



IJiMifii) ® \J 2 M 2 n 2 ) 



(21) 



The dressed state (DS) van der Waals interaction 
i?Pf (R) is calculated from the matrix elements of Eq. (|20|) 
in the dressed state basis, 
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TABLE IV. XY (X,Y=Na,K,Rb,Cs) calculated CCSD(T) electrostatic and TD-DFT dispersion+induction van der Waals 
coefficients, W^^ L M , for unique combinations of L1L2AI. All values are presented in atomic units and calculated at the 
equilibrium bond length r e listed in Table U and [n] denotes x!0 n . 



n L\ Li M 


NaK 


NaRb 


NaCs 


KRb 


KCs 


RbCs 


Electrostatic: 














3110 


-2.470[0] 


-3.651[0] 


-6.813[0] 


-0.125[0] 


-1.164[0] 


-0.432 [0] 


3111 


1.235[0] 


1.826[0] 


3.406[0] 


0.063[0] 


0.582[0] 


0.216[0] 


4210 


-3.528[1] 


-2.818[1] 


-1.377[1] 


-1.134[1] 


-2.960[1] 


-2.221 [1] 


4211 


1.176 [1] 


9.393[0] 


4.590[0] 


3.785[0] 


9.867[0] 


7.403[0] 


5220 


6.717[2] 


2.900[2] 


3.710[1] 


1.375 [3] 


1.004 [3] 


1.524 [3] 


5221 


-1.493 [2] 


-6.444[1] 


-8.245 [0] 


-3.056[2] 


-2.231 [2] 


-3.386[2] 


5222 


9.329[0] 




U.010[UJ 


1 Qinfii 




Z. 11D[1J 


5310 


1.175 [2] 


7.346[1] 


5.962[1] 


6.909[1] 


1.955[2] 


1.106[1] 


5311 


-2.937[1] 


-1.836[1] 


-1.490[1] 


-1.727[1] 


-4.887[1] 


-2.764[0] 


Dispersion+induction: W^J L M 












6000 


7.777[3] 


8. 680 [3] 


1.233 [4] 


1.354[4] 


1.726 [4] 


1.921 [4] 


6200 


5.519[2] 


7.837[2] 


3.327[2] 


2.001 [3] 


2.375[3] 


2.909[3] 


6220 


9.762[2] 


1.223[3] 


2. 694 [3] 


1.028[3] 


1.826[3] 


1.857[3] 


6221 


-2.169 [2] 


-2.717[2] 


-5.986 [2] 


-2. 284 [2] 


-4.059[2] 


-4.127[2] 


6222 


2.712[1] 


3.397[1] 


7.482[1] 


2.855[1] 


5.073[1] 


5.159[1] 


7100 


1.157[4] 


2. 268 [3] 


-2.114[4] 


6.069[3] 


1.457 [4] 


1.755 [4] 


7210 


-1.657[1] 


-4.434[3] 


-2.312[4] 


9.929[2] 


-8.871 [2] 


2. 795 [3] 


7211 


2.762[0] 


7.391 [2] 


3.853[3] 


-1.655 [2] 


1.479 [2] 


-4.658[2] 


7300 


9.136 [3] 


4.795[3] 


-1.280 [3] 


4.176[3] 


1.234 [4] 


1.231 [4] 


7320 


4.931 [3] 


1.142 [3] 


-9.091 [3] 


2. 772 [3] 


7.134[3] 


8.669[3] 


7321 


-8.219[2] 


-1.903[2] 


1.515[3] 


-4.621 [2] 


-1.189 [3] 


-1.445 [3] 


7322 


5.870[1] 


1.360[1] 


-1.082 [2] 


3.300[1] 


8.492[1] 


1.032[2] 


8000 


1.444 [6] 


1.928[6] 


3.016[6] 


3. 734 [6] 


5.391 [6] 


5.667[6] 


8200 


8.920[5] 


1.503 [6] 


2. 766 [6] 


2. 606 [6] 


4.295[6] 


3.923[6] 


8220 


3.296[5] 


6.003[5] 


1.466[6] 


7.941 [5] 


1.494 [6] 


1.296[6] 


8221 


-4.768 [4] 


-8.470[4] 


-2.042[5] 


-1.162 [5] 


-2.138[5] 


-1.878[5] 


8222 


3. 368 [3] 


4.596[3] 


9. 349 [3] 


9.098[3] 


1.356 [4] 


1.351 [4] 


8400 


6. 476 [4] 


6.275[4] 


1.107[5] 


2.931[5] 


3.258[5] 


4.238[5] 


8420 


5.979[4] 


7.214[4] 


1.962 [5] 


2.200[5] 


2.678[5] 


3.389[5] 


8421 


-8.081[3] 


-9.702 [3] 


-2.625[4] 


-2.970[4] 


-3.612[4] 


-4.566[4] 


8422 


4.380[2] 


5.183[2] 


1.380 [3] 


1.603 [3] 


1.945 [3] 


2. 453 [3] 



Ji.Mi J2M2 Li,m Ll 
J[,M[ J 2 ,M' 2 L 2 ,m L2 



(^j^-a^^^^^-ft^^A^^^ [(2Ji + 1)(2J , + 1)(2J2 + 1)(2J / + 1)(2Li + 1)(2ia + 

,( J'\ L x Ji \ / J[ L X Ji \ / J' 2 L 2 J-2 \ fj-2 L 2 J 2 \ ST^ ( L l L 2 L \ V I f}\ pLF ( 

\M[ m Ll -MJ Ifii -nj \M 2 m L2 -M 2 Ul 2 -fi 2 J 2^ \m Li m La m J 1 ^K^l^lK 

L.7UL 



/2 

(R), 



(22) 



where 

Pm%i( f ) = a 'U F ) a fu( F ) (23) 
is the coupled rotational state density of molecule i and 



YL m (f) is a spherical harmonic [38j]. In addition to the 
transformation of the van der Waals interaction energy 
as given by Eq. (j2"2"|) . it is useful to have the dressed static 
moment, (Qf s ), of a given molecule. For molecule i, this 



7 



a-e CCSD(T)-fl2a/QZVPP 

— vdW TD-DFT PBEO 

-- C 3 (9)/R 3 +C 5 (8)/R 5 - C 6 ' s °/R 6 



Collinear 




40 50 
R (a.u.) 

FIG. 2. Comparison between the TD-DFT van der Waals 
surface, evaluated at both collinear (6 = 0) and aligned (0 = 
7r/2) geometries, and both the electrostatic + isotropic C§ 
and a fully ab initio curve for LiNa+LiNa. 



is readily obtained to be 

(Qm S ) = (JlMln'iliOilQioMJiMiQi) = 

E ((2J< + i)(2J| + 1)) 1/2 



X(-1) 



m; o -M 



fii -Q, 



(Q/o>- (24) 



Low-Field Solution 



In the low-field limit, coupling between rotational 
states can be limited to just two states, allowing Eq.© 
to be solved analytically (note that = and M = 0, 
as discussed below in Sec. IVI[) . ^From this it is possible 
to obtain general expressions for the expectation value of 
the static and alignment moments as a function of the ap- 
plied field. Transforming to the unitless field parameter 
£ = £ F, with Co = V/2B, the low-field limit is defined 
by £ < 1. With this approximation, the dressed state 
dipolc and quadrupolc moments can be shown to be 



(Q?o s )(0 = (Qio)^^ 



and 



(QgfXO = {Q70) 



2e 



15(1 + £ 2 ) 



(25) 



(26) 



respectively, while the octupole moment has no two state 
contribution by symmetry. The orientation moment 
(cos 9) is given trivially by 



(cos0>(O = (Q$)(0/{Qw), 



(27) 



while alignment (cos 2 9) can be calculated by noting that 
cos 2 # = !(l+2Ci i o(#)) (where C^ m {fi) is a Racah spher- 
ical harmonic), providing the expression 



(cos' 



15 + 19<f 

45(1 +ey 



(28) 



So long as the number of coupled states is dominated 
by the first two states and £ < 1, these approximate 
formula are accurate to a few percent. In Table [V] we 
have evaluated £o for all the heteronuclear alkali diatoms 
from the spectroscopic data in Table HI It is also possible 
to evaluate Eqs. (|2"2"|) and ([2"B")) using the two state low- 
field approximation. Following the prescribed method 
discussed above, the low field-dressed-state van der Waals 
potential can be written to leading order as 



(i) 



R 3 

(2) 
6000 



R 5 



(2) 



8000 



R e 



R 8 



(29) 



Here 9 F is the angle between the inter-molecular vec- 
tor R and the field vector as illustrated in Fig. [T^b). 
The dipolc-dipole and quadrupolc-quadrupolc contribu- 
tions are (up to order £ 5 ) 



■320\"F,C) - (QlO 

x(i-; 



3v^ + 6£ 2 + 4\/3£ 3 + 4£ 4 



27(1 +ey 2 

f), 



(30) 



and 



75(1 +e) 2 

x (3 -30 cos 2 9 F + 35 cos 4 9 F ) , (31) 



respectively (note that there is no dipole-octupole con- 

(2) 

tribution in the two state approximation) while !F„ oo i s 
the isotropic dispersion+induction coefficient (see Tables 
IIIII and lIVp . The anisotropic terms contribute less than 
a percent to the interaction energy and can be safely ne- 
glected. 



V. ELECTRONIC STRUCTURE 
CALCULATIONS 

The ab initio calculation of van der Waals coefficients, 
and more generally multipole polarizabilities, requires 
special care in both the basis set and level of theory used 
[4ll - |46| . Electrostatic moments similarly require careful 
consideration of the theoretical method, though the ba- 
sis set dependence is less severe [13] . For all calculations 
in this work we use the Karlsruhe def2-QZVPP [48| basis 
set augmented with two additional even tempered diffuse 
spdf functions designed to accurately describe higher or- 
der static polarizabilities [3l[ . The Karlshruhe def2 basis 
sets are available for nearly the entire Periodic Tabic, 
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FIG. 3. Dressed state electrostatic moments, (Qth }> °f various heteronuclear alkali diatomic molecules as a function of an 
external DC electric field. 



and are known for both their robustness and good cost- 
to-performance ratio in large molecular Hartree-Fock and 
density functional theory calculations. As such they re- 
main attractive for use in calculations that involve many 
different atoms across the period table. 

As has been demonstrated previously, the use of time 
dependent density functional theory [43, E(| is a cost ef- 
fective and accurate way to calculate multipole transition 
moments and excitation energies for diatomic molecules. 
We chose to limit our calculations in this work to only 
include the PBEO functional for simplicity, however for 
various cases it was observed that the B3PW91 func- 
tional also provides consistent results. The electrostatic 
moments were calculated using coupled cluster theory 
including all singles, doubles and perturbative triples 
(CCSD(T)) [U using a two step finite field method (with 
field spacings of 10~ 6 a.u.). Core- valence and core-core 
correlation energy was accounted for by including the 
inner valence s and p electrons in the CCSD(T) calcula- 
tions, while for the TD-DFT computations all electrons 
not replaced by an ECP are implicitly correlated. All 



TD-DFT calculations were done using a locally modi- 
fied version of the GAMESS [HI, suite of programs; 
the CCSD(T) finite field calculations were done using the 
MOLPRO [54j quantum chemistry program package. For 
further details on the methodology used in evaluating the 
transition dipole moments and excitation energy we re- 
fer to our previous paper on homonuclear alkali diatomic 
molecules f3lj . 



VI. COMPUTATIONAL RESULTS AND 
DISCUSSION 

The leading order term of the long range expansion 
Eq.@, and thus the longest ranged interaction in the 
series, involve products of the electrostatic moments of 
each monomer, and for dipolar molecules, it is the dipole- 
dipole R~ 3 term. The dipole-dipole scatte ring [f35[ and 
applications of dipole-dipole interactions [l4[ are well 
studied in the literature, however higher order terms can 
be necessary for accurately describing intermediate inter- 
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FIG. 4. DC field coupled van der Waals curves (Eq. (J22J ) of 40 K 87 Rb for both low and high fields as well as the approximate 
two-state van der Waals curve (Eq. (J29J ) . Here F ~ 20 kV/cm is the intermediate field strength where more than two rotational 
states begin to strongly couple. 



molecular distances [3l( and are often neglected if only for 
a lack of available data. Inclusion of just the quadrupolc- 
quadrupole interaction to a dipole-dipole model can in- 
troduce significant changes in the form of potential en- 
ergy barriers for collinear geometries (9\ = 6 2 = <f> = 0) 
at long range (56|. It is possible to estimate whether the 
inclusion of higher order electrostatic terms could lead 
to a barrier by introducing the outer zero energy turning 
point, R q , which occurs when the R~ 5 repulsion over- 
comes the attractive R~ 3 dipole-dipole force. Keeping 
only the two leading terms in Eq.(|9]) and setting 



E int (R q , 0,0,0) 
we obtain 



Rq = 



(Qw) 2 , 3(Q 2 o) 2 -4(Q 10 )(Q 3 o) 



R 5 g 



V / 3(Q 2 o) 2 -4(Q 10 )(Q 



:;u/ 



(Qi 



= 0, 
(32) 

(33) 



When this outer turning point is sufficiently long range 
(R q > 20 a.u.) the introduction of these higher order 
terms can be important and lead to long range barri- 
ers [56|, and thus should be examined in further detail. 
As such, we have calculated the ab initio electrostatic 
dipole, quadrupole and octupole moments (higher order 
moments do not contribute up to i?~ 5 in the long range 



expansion: see Sec. |V] for details on the methodology 
used). In Table U we present our calculated static mo- 
ments, the outer turning point R q for each system, as 
well as various dipole and quadrupole moments found 
in the literature. Our computed static dipole moments 
agree closely with both the valence full configuration 
interaction results of Aymar et al. [34[ across all the 
molecules investigated, and the CCSDT (CCSD with all 
triples) results of Quemener et al. Q for the highly polar 
LiX (X=Na,K,Rb,Cs) species. Other than the CCSD(T) 
quadrupole moment of Zemke et al. [27j (with which we 
compare well), little to no published quadrupole values 
exist for the heteronuclear alkali diatoms. It has been 
demonstrated for the homonuclear alkali diatoms that 
the finite field CCSD(T) higher order static moments 
compare well with other methods (3ll . [47j; similar ac- 
curacy is anticipated for the heteronuclear species. 

Dispersion and induction contributions to the van der 
Waals series are proportional to products of the dipole, 
quadrupole and octupole polarizabilities. As such we 
have calculated and presented in Table [TT| the dipole 
and quadrupole static polarizabilities with comparisons 
to some of the existing literature (octupole static polar- 
izabilities are not listed, but are available upon request). 
As discussed previously [3l|, the n-aug-def2-QZVPP ba- 
sis sets are well converged for computation of static polar- 
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izabilitics of homonuclcar alkali diatoms up to octupolc 
order, and we find the same is true for the heteronu- 
clear species. Density functional methods are known to 
provide average static polarizabilities to within five to 
ten percent of experimental or highly correlated results 
[57L [58^ . Furthermore some variance is expected in the 
parallel (ano) polarizability as all computations are done 
at the experimental (or theoretical where necessary) equi- 
librium bond length, and it is well known that the polar- 
izability is sensitive to the internuclear separation in the 
alkali diatoms [33j . It is expected that the perpendicular 
polarizability (am) should agree much more closely with 
other methods, which we find to be the case as illustrated 
in Table IU 

Van der Waals dispersion and induction coefficients of 
the heteronuclear alkali diatoms are sparsely given in 
the literature. Currently only a few values exist and 
are restricted to isotropic contributions (corresponding 

to M^gooo )■ The only systematically calculations are for 
the LiX species Q. In Table [TT| we note the reason- 
able agreement between our reported TD-DFT isotropic 



C 6 = W, 



2.DIS 
60oq {- 



dispersion coefficients and the Tang-Slater- 



Kirkwood [36| values from Qucmcncr et al. [8[ for the 
LiX species. Additionally Kotochigova [32[ has calcu- 
lated, using multi-reference configuration interaction the- 
ory, the isotropic and anisotropic dispersion coefficients 
of order R~ 6 for both KRb and RbCs. However, these 
values contain non-Born-Oppcnheimer contributions and 
so arc not directly comparable to our numbers; be- 
cause of this we have not included these values in Ta- 
ble [TTJ To determine the accuracy of the van der Waals 
coefficients calculated here, we have computed ab ini- 
tio curves for LiNa+LiNa at two different geometries 
using the CCSD(T)-F12a/QZVPP (explicitly correlated 
CCSD(T)) level of theory j53.[60|. These ab initio curves 
are plotted in Fig. [2] along with the electrostatic plus 
isotropic dispersion approximation and the the van der 
Waals curves of this work including all anisotropic terms 
through R~ 8 . As can be seen, the TD-DFT van der 
Waals curves agree to a few cm -1 with the ab initio re- 
sults [6l[, while the isotropic curves fail completely in 
the intermediate-range (it should be noted that for the 
collinear case of LiNa+LiNa the isotropic curves do not 
turn over at all and predict an infinite repulsive wall). In 



Tablcs[ITT]and|TV]wc have listed the W. 



(1,2) 



coefficients 

for all of the heteronuclear alkali diatoms, including all 
terms up through order i?~ 8 . 

In evaluating the field coupling and alignment of the 
various alkali diatomic molecules, the rotational wave- 
function expansion is greatly simplified by making use 
of the initial premise that the molecules are in the ro- 
vibrational ground state. As such S7 = and M = 
(the use of a DC external field will not mix different 
M values), reducing both Eqs. (f2"2"]) and (j2~4"]l significantly. 
In Fig. [3] we have plotted the DC field dressed elec- 
trostatic moments as a function of the external field 
strength. While the very high field strengths in Fig. [3] 
are generally experimentally challenging, it is illustrative 



TABLE V. Tabulated values of the field strength coefficient 
£0 = T>/2B using the spectroscopic and electrostatic constants 
from Table IT] All units are in cm 2 /kV. 





23 Na 


39 K 


95 Rb 


133 Cs 


'Li 


0.0114 


0.116 


0.159 


0.246 


23 Na 




0.253 


0.413 


0.684 








0.141 


0.529 


95 Rb 








0.635 



to show how difficult it is to obtain both strong orienta- 
tion ((cos#) > 0.85) and alignment ((cos 2 #) > 0.85) in 
molecules with small rotational constants, regardless of 
the strength of the dipole moment. It is also instructive 
to examine the low-field strengths of Fig. [3l where the 
linear trend of each curve on the log-log scale shows the 
general scaling of the static moments as a function of the 
external field as discussed in Sec. IIVBI In Fig 0] we have 
evaluated Eq. ([22|) for KRb (KRb is chosen for its medium 
strength dipole moment and large rotational constant) at 
various DC field strengths. The difference between the 
low and high field strengths is easily identified by the 
change in behavior from most similar to the field free 
case (e.g. isotropic contributions dominate the interac- 
tion potential) to the regime where the dressed state van 
der Waals interaction energy more closely resembles the 
molecule-fixed frame van der Waals potential (e.g. when 
electrostatic contributions become key). This high field 
strength regime is more quantitatively defined when both 
(cos 6) and (cos 2 9) is greater than 0.9 (which corresponds 
to roughly 7 strongly coupled rotational states). Also in 
FigU] the approximate two-state model of Eq. ([2"9")l can 
be seen to agree very well with the fully coupled equa- 
tions in the low- field limit. Fully field-coupled potentials 
for the other heteronuclear molecules listed in this work 
have been calculated, and are available upon request. 



VII. CONCLUSIONS 

This work completes our systematic TD-DFT compu- 
tation of the alkali diatomic species by computing ac- 
curate multipole electrostatic moments and anisotropic 
van der Waals coefficients for the heteronuclear alkali 
diatomic species. The multipole electrostatic moments 
were computed using a finite field treatment of the 
CCSD(T) molecular energy employing the augmented 
Karlsruhe dcf2-QZVPP basis set and found to produce 
excellent agreement with the existing literature. Excita- 
tion energies and multipole transition moments were cal- 
culated using TD-DFT and the same augmented QZVPP 
basis set. Static polarizabilities as well as van der Waals 
induction and dispersion coefficients were evaluated using 
the sum over states approach and found to be consistent 
with the existing literature. Using the simple form of 
Eq. © and the values from Tables Mil and IIV1 it is possi- 
ble to completely characterize the long range interaction 
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between two heteronuclear alkali diatoms up through or- 
der R~ 8 . A sample FORTRAN program for evaluating 
Eq.(E3) is included in the supplemental material of Rcf. 
[3 1| or upon request to the authors. 

The transformation of the van der Waals series for lin- 
ear molecules from the molecule-fixed frame to the lab- 
fixed frame was described. This was followed with the 
computation of the dressed state electrostatic moments 
as a function of an external DC electric field. It was noted 
that in the low field limit, the coupling of the molecule to 
the external field can be approximated by only consider- 
ing two rotational states. With this in mind, the orien- 
tation and alignment of the molecule as a function of the 
applied field can be approximated using only molecular 
spectroscopic constants by Eqs. (p3l) - (|28|) . which are valid 
for field values F < 2B/D. We have also illustrated the 
effects of an external DC electric field on the intermolec- 



ular potential by evaluating Eq. (|22"j) for 40 K 87 Rb at a va- 
riety of field strengths. It can be seen then that introduc- 
ing rotational state coupling leads to a richer interaction 
phase space beyond the usual isotropic approximations. 
Finally a two-state approximation of the dressed-state 
long range potential (see Eq. (|22]) ) has been derived and 
given by Eq. (|29|) in terms of molecular spectroscopic con- 
stants and isotropic van der Waals coefficients 
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